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Vega  is  a  rapidly  rotating  star 

D.  M.  Peterson1,  C.  A.  Hummel2’3,  T.  A.  Pauls4,  J.  T.  Armstrong4,  J.  A.  Benson",  G.  C.  Gilbreath4,  R.  B.  Hindsley4, 
D.  J.  Hutter5,  K.  J.  Johnston",  D.  Mozurkewich6  &  H.  R.  Schmitt4,7 


Vega,  the  second  brightest  star  in  the  northern  hemisphere,  serves 
as  a  primary  spectral  type  standard1.  Although  its  spectrum  is 
dominated  by  broad  hydrogen  lines,  the  narrower  lines  of  the 
heavy  elements  suggested  slow  to  moderate  rotation,  giving 
confidence  that  the  ground-based  calibration  of  its  visible  spec¬ 
trum  could  be  safely  extrapolated  into  the  ultraviolet  and  near- 
infrared  (through  atmosphere  models2),  where  it  also  serves  as  the 
primary  photometric  calibrator.  But  there  have  been  problems: 
the  star  is  too  bright  compared  to  its  peers3  and  it  has  unusually 
shaped  absorption  line  profiles,  leading  some4,5  to  suggest  that  it  is 
a  distorted,  rapidly  rotating  star  seen  pole-on.  Here  we  report 
optical  interferometric  observations  that  show  that  Vega  has  the 
asymmetric  brightness  distribution  of  the  bright,  slightly  offset 
polar  axis  of  a  star  rotating  at  93  per  cent  of  its  breakup  speed.  In 
addition  to  explaining  the  unusual  brightness  and  line  shape 
peculiarities,  this  result  leads  to  the  prediction  of  an  excess  of 
near-infrared  emission  compared  to  the  visible,  in  agreement  with 
observations6,7.  The  large  temperature  differences  predicted 
across  its  surface  call  into  question  composition  determinations, 
adding  uncertainty  to  Vega’s  age  and  opening  the  possibility  that 
its  debris  disk8  could  be  substantially  older  than  previously 
thought9,10. 

Single  baseline  Michelson  stellar  interferometers  measure  complex 
‘visibilities’11,  usually  recorded  as  amplitudes  and  phases,  which  are 
related  to  the  intensity  distribution  of  the  target  through  a  Fourier 
transform.  Even  though  the  phases  have  been  shown  to  be  very 
sensitive  to  asymmetries  in  the  intensity  distribution,  they  are  badly 
corrupted  by  the  atmosphere  and  have  been  little  used  in  the  optical. 
But  for  closure  phases,  the  data  we  focus  on  here — which  are 
obtained  by  summing  the  phases  measured  on  each  baseline  of  a 
triangle  in  an  interferometric  array  such  as  the  Navy  Prototype 
Optical  Interferometer12  (NPOI) — the  atmospheric  contribution 
cancels.  The  use  of  closure  phase  in  the  radio13  has  enabled  a  dramatic 
gain  in  dynamic  range  of  interferometric  images  made  from  multi¬ 
antenna  arrays.  In  the  optical12,14,  where  the  phase  errors  can  reach 
100  waves  on  long  baselines,  the  technique  enables  the  use  of  phase 
information  in  any  guise.  As  the  observations  reported  here  were 
made  just  with  a  three-telescope  array,  the  application  of  imaging 
techniques  was  not  justified  and  we  relied  instead  on  fitting  Roche 
models  to  the  closure  phase  data. 

The  application  of  Roche  spheroids  to  rotating  stars  was  worked 
out  80  years  ago15.  Assuming  solid  body  rotation  and  a  point  mass 
gravitational  potential,  a  rotating  star  will  adopt  the  figure  of  a  Roche 
spheroid.  Conservation  of  energy  through  surfaces  of  constant 
potential  leads  to  the  prediction  that  when  the  energy  is  transported 
by  radiation  the  amount  transported  will  vary  over  the  surface  in 
proportion  to  the  effective  gravity16  (that  is,  gravity  minus  local 
centrifugal  terms).  Near  breakup,  the  effective  gravity  near  the 
equator  can  become  quite  small,  leading  to  the  prediction  of  a 


large  drop  in  the  local  temperature  with  a  corresponding  decrease 
in  brightness,  an  effect  referred  to  as  ‘gravity  darkening’.  Rapidly 
rotating  stars  seen  at  intermediate  inclinations  are  therefore  expected 
to  display  asymmetric  intensity  distributions.  Altair17,18  proved  to  be 
the  first  major  test  of  this  theory  in  an  isolated  rotating  star,  where  the 
theory  succeeded  to  a  high  degree  in  describing  a  very  non-trivial 
brightness  distribution. 

The  observations  considered  here  were  obtained  during  late  May 
and  early  June  2001  on  the  same  nights  as  those  of  Altair  previously 
reported17,19,  where  Vega  served  as  a  check  star  (the  Vega  data  in 
machine  readable  form  are  in  a  separate  file  in  the  Supplementary 
Information).  The  observations  and  much  of  the  data  reductions 
are  exactly  as  described  for  Altair17,  which  should  be  consulted 
for  details.  Issues  specific  to  the  Vega  data  are  described  in  the 


Table  1  |  Vega  model  and  derived  parameters 


Quantity 

Value 

Errors  (s.d.)* 

co  —  (2/  (2  bT 

0.926 

±0.021 

Op  (mas)t 

2.767 

0.037 

T p  (K)t 

9,988 

61 

i  (deg.)t 

4.54 

0.33 

PA  (deg.)t 

8.6 

2.7 

veq  (km  s_1) 

274 

14 

^ eq,B-i-  (km  S-1) 

356.1 

2.4 

0  (d_1) 

1.884 

0.081 

ObX  (d_1) 

2.034 

0.041 

Te q  (K) 

7,557 

261 

R  p(Rq) 

2.306 

0.031 

R eq  (Rq) 

2.873 

0.026 

0min§  (mas) 

3.441 

0.031 

0ma*§  (mas) 

3.446 

0.031 

log  1(1.0) 

1.544 

0.018 

log  SP  (cm2s  2) 

4.074 

0.012 

log  Seq  (cm2S~2) 

3.589 

0.056 

M||  (M0) 

2.303 

0.024 

Feffll  (K) 

9,306 

86 

Age ||  (Myr) 

386 

16 

Six  quantities  are  needed  to  uniquely  define  the  Roche  model  of  a  star17:  the  ratio  of  the 
angular  rotation  to  that  of  breakup,  co  =  (2/(2 b,  the  inclination  (or  tilt)  of  the  rotational  axis  to 
the  line  of  sight,  /,  the  position  angle,  PA,  of  the  pole  on  the  sky,  the  radius  of  the  polar  axis, 

R p  (or  equivalently,  using  the  parallax,  the  polar  angular  diameter,  0P),  the  effective 
temperature  at  the  pole,  7P,  and  the  surface  gravity  at  the  pole,  gp,  or  equivalently,  the  mass. 
It  is  then  possible  to  calculate  the  radius,  R(</>),  of  the  star  for  a  given  stellar  latitude,  </>,  the 
effective  gravity  and  the  local  temperature  (7(0)).  For  the  spectral  calculations  we  have 
adopted  the  ATLAS  model  atmospheres28  and  in  particular  the  Van  Hamme29  limb-darkening 
parameterization  of  that  grid.  Other  parameters  include  linear  rotational  velocities,  v,  the 
luminosity,  L,  surface  gravities,  g,  the  mass,  M,  and  effective  temperature,  7eff,  the  last 
referring  to  the  entire  non-rotating  star.  Subscripts  eq  and  p  specify  quantities  evaluated  at 
the  equator  and  pole,  respectively. 

‘Uncertainties  due  to  the  parallax  have  not  been  included  in  the  errors. 

fThe  parameters  derived  from  the  model  fit.  A  mass  of  2.30  M©  was  assumed  in  the  fit. 

:j:  Rotating  at  breakup  but  with  the  same  mass  and  polar  radius. 

§0min,max  are  the  minimum  and  maximum  projected  angular  diameters. 

||  The  parameters  of  a  non-rotating  star  from  the  Padova  grid26  that  would  reproduce  the 
(corrected)  luminosity  and  polar  radius17. 
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Table  2  |  Predicted  photometric  'excesses'  for  Vega 


Band 

Wavelength  (nm) 

Excess  (mag) 

R 

0.67 

-0.016 

1 

0.86 

-0.029 

J 

1.12 

-0.052 

K 

2.14 

-0.072 

L 

3.69 

-0.072 

The  differential  excesses  for  the  model  of  Vega  compared  to  a  non-rotating  model  (both 
solar  composition)  which  matches  the  V  magnitude  with  an  angular  diameter  of  3.24  mas 
(ref.  27).  Until  the  question  of  composition  is  resolved,  these  values  should  be  considered 
indicative. 


Supplementary  Information.  We  focus  here  on  the  closure  phase  data 
taken  by  NPOI  on  25  May  2001,  as  they  are  the  most  extensive  and  the 
highest  quality  data  of  that  run.  We  augment  these  with  the  V-band 
magnitude2,  V  =  0.026  ±  0.008,  as  an  additional  observable.  The 
model  was  fitted  by  enforcing  the  usual  minimum  x 2  metric  using 
the  Levenberg-Marquardt  algorithm20. 

The  parameters  from  the  initial  reduction  are  given  in  column  2  of 
Supplementary  Table  SI.  The  projected  (at  inclination  i )  equatorial 
velocity,  veq,  for  this  model  was  predicted  to  be  veq  sin  i  ~  lSkms^1, 
a  bit  below  the  —21.8  ±0.1  km  s  1  found  from  detailed  profile  fits 
using  a  rotating  model21.  Although  it  is  not  so  far  off  given  how  close 
the  star  is  to  pole-on,  we  feel  the  projected  velocity  is  sufficiently  well 


known  that  it  should  also  be  incorporated  in  the  fit.  We  have 
therefore  added  as  an  ‘observable’  veqsinz  =  22  ±  1.0  km  s  1  (see 
Supplementary  Information)  and  refit  the  data  (columns  3  and  4  of 
Supplementary  Table  SI). 

The  Roche  model  provides  a  good  fit  to  the  augmented  data  set. 
The  model  parameters  fitted,  a  number  of  derived  quantities,  and  an 
estimate  of  what  the  main  parameters  of  the  star  Vega  would  be,  were 
it  not  rotating,  are  given  in  Table  1 .  The  quality  of  the  fit  is  illustrated 
in  Fig.  1,  and  a  false-colour  rendering  of  the  Vega  model  is  shown  in 
Fig.  2.  As  Vega  has  so  many  roles  in  astronomy,  the  ramifications  of 
this  result  are  extensive.  We  summarize  some  of  the  most  important 
below. 

Rotation  can  affect  the  gross  spectral  distribution  of  a  star,  an  issue 
of  considerable  import  given  Vega’s  role  as  the  primary  flux  calibrator 
in  the  ultraviolet,  visible  and  near-infrared.  To  estimate  the  size  of 
these  effects,  we  have  calculated  the  changes  in  the  fluxes  from  the 
rotating  model  compared  to  the  static  case  that  would  be  measured 
through  a  series  ofbroadband  filters  (Table  2).  As  can  be  seen,  there  is 
a  significant,  systematic  increase  in  the  infrared  emission  from  the 
rotating  model.  This  is  understood  as  being  due  to  the  large  amount 
of  surface  area  predicted  to  be  at  relatively  low  temperatures. 

There  is  an  extensive  literature  on  the  possibility  and  extent  of  an 
‘infrared  excess’  in  the  Vega  spectrum6,7.  As  the  issue  is  critical  to  so 
much  of  astronomy,  sides  have  been  strongly  taken.  Observations  of 
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Figure  1  Roche  model  fits  to  the  closure  phase  data  taken  on  25  May 
2001.  The  observations  (open  circles),  estimated  errors  (bars,  standard 
deviations)  and  model  calculations  (solid  lines)  are  shown  for  each  scan 
(labelled  by  hour  angle,  HA).  Residuals  (observed,  O,  minus  calculated,  C) 
are  shown  below  each  of  the  scans  for  clarity.  The  phases  for  the  individual 
baselines,  o„  and  thus  the  closure  phases,  take  on  only  two  values,  0°  or  180°, 


if  an  object  is  centro-symmetric.  Closure  phase  measurements  showing 
departures  from  this  simple  ‘abrupt  transition’  behaviour  provide 
potentially  very  sensitive  measurements  of  asymmetry  in  an  object.  The  soft 
transition  at  the  points  of  the  180°  phase  changes  here  give  a  clear  signal  of 
the  asymmetry  in  the  intensity  distribution.  (Note  the  scale  change  for  the 
residuals  of  the  last  scan.) 
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North 


Figure  2  |  A  false-colour  model  of  Vega  as  it  appears  from  Earth.  (Blue  is 
bright,  red  is  faint,  and  the  orange  '+’  is  the  subsolar  point.)  The 
temperature  drops  more  than  2,400  K  from  pole  to  equator,  creating  an  18X 
drop  in  intensity  at  500  nm.  Limb-darkening  in  a  non-rotating  model 
predicts  only  a  5-fold  drop  in  intensity.  Although  the  projected  outline  is 
almost  perfectly  circular,  the  polar  diameter  is  only  80%  of  the  equator. 
Dec.,  declination;  RA,  right  ascension. 

other  A  stars  showed  that  Vega’s  colours  appear  sensibly  normal7, 
which  led  some  authors7,22  to  wonder  whether  the  problem  was  with 
the  model  atmospheres.  Others23  argued  that,  because  the  hydrogen 
absorption  coefficient  completely  dominates  the  spectra  of  A  stars 
and  is  so  well  known,  it  was  unlikely  that  the  model  atmospheres  were 
wrong,  an  argument  that  seems  to  have  carried  the  day.  Rapid 
rotation  provides  a  simple  resolution  to  this  controversy — Vega  is 
best  modelled  as  a  composite  of  model  atmospheres.  The  star  should 
have  ‘normal’  infrared  colours,  as  most  normal  A  stars  are  rapid 
rotators,  and  there  is  no  need  to  question  the  quality  of  the  individual 
atmosphere  models. 

One  important  aspect  to  this  model  is  the  near  pole-on  orientation 
of  the  rotational  axis.  Vega  is  surrounded  by  a  large  infrared- emitting 
disk  of  material,  a  ‘debris  disk’,  which  presents  an  essentially  circular 
profile.  We  do  not  expect  perfect  coupling  between  orbiting  material 
and  the  central  star.  But,  if  the  poles  of  the  disk  and  Vega  coincide  and 
the  disk  is  thin,  we  would  predict  0.3%  flattening,  which  is  unlikely  to 
be  detectable,  as  seems  to  be  the  case24. 

Related  to  the  orientation  of  the  pole  is  the  inferred  equatorial 
velocity,  veq  =  272  km  s  1  (Table  1).  Among  the  Vega-like  stars,  Vega 
itself  has  been  anomalous  in  displaying  a  very  low  projected 
rotational  velocity10.  The  present  determination  strongly  supports 
the  view  that  the  Vega-like  stars  are  rapid  rotators,  consistent  with  the 
large  amounts  of  angular  momentum  in  the  surrounding  dust 
clouds. 

Vega’s  rotational  state  also  affects  inferences  about  its  debris  disk.  It 
is  well  known  that  rotation  results  in  the  apparent  brightness,  and  in 
turn  the  deduced  luminosity,  being  inclination-dependent25.  Song 
et  al. 10,  for  example,  have  gone  to  some  lengths  to  characterize  the 
uncertainty  introduced  by  this  effect  in  the  Vega-like  stars.  It  is 
straightforward  to  calculate  the  total  luminosity  of  a  Roche  spheroid 
(Table  1).  Further,  one  can  apply  small  corrections17  to  the  lumi¬ 
nosity  and  polar  radius  to  obtain  the  corresponding  values  that 
would  apply  to  a  non-rotating  star  of  the  same  mass.  Using  the 
Padova26  models,  we  derive  a  mass  M  =  2.303  Mq  (where  Mq  is  the 
solar  mass)  and  an  age  of  386  ±  16Myr,  on  the  high  side  of  recent 
estimates  of  354±g°Myr  (ref.  10)  and  347l33Myr  (ref.  9). 

Unfortunately,  Vega’s  composition  enters  the  age  determination 
rather  critically.  The  star  is  currently  viewed  as  underabundant  in 
heavy  elements  compared  to  the  Sun27,  [Fe/H]  ~  —  0.5,  which 
clearly  needs  to  be  revisited,  given  the  2,400  K  temperature  drop 


now  predicted  across  its  surface.  Previous  age  determinations  have 
implicitly  or  explicitly  assumed  solar  metallicity,  the  argument  being 
that  compositional  peculiarities  among  the  slowly  rotating  stars  in 
this  part  of  the  Hertzsprung-Russell  diagram  are  probably  limited  to 
the  outer  envelope,  and  normal  composition  models  are  therefore 
appropriate  for  estimating  bulk  properties.  In  our  interpretation  that 
argument  fails,  because  at  these  rotational  velocities  meridional 
circulation  will  keep  the  bulk  of  the  star  well  mixed. 

Proceeding  as  above,  we  estimate  M=  2.11  Mq  and  an  age  of 
572 Myr,  using  the  metallicity  Z  =  0.008  ([Fe/H]  ~  —0.4)  Padova 
models26.  It  is  clear  that  a  full  abundance  analysis  incorporating 
rotation  needs  to  be  performed  to  remove  this  uncertainty.  In  the 
meantime  this  range,  386-572  Myr,  is  probably  a  more  realistic 
estimate  of  the  uncertainty  in  the  evolutionary  age  of  Vega  and  its 
debris  disk. 
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